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ABSTRACT 

We present sensitive radio continuum imaging at 1.4 GHz and 4.9 GHz of seven high redshift 
QSOs selected for having a 240 GHz continuum detection, which is thought to be thermal dust emis- 
sion. We detect radio continuum emission from four of the sources: BRI 0952—0115, BR 1202—0725, 
LBQS 1230+1627B, and BRI 1335-0417. The radio source in BR 1202-0725 is resolved into two com- 
ponents, coincident with the double mm and CO sources. We compare the results at 1.4 GHz and 240 
GHz to empirical and semi-analytic spectral models based on star forming galaxies at low redshift. The 
radio-to-submm spectral energy distribution for BR 1202-0725, LBQS 1230+1627B, and BRI 1335-0417 
are consistent with that expected for a massive starburst galaxy, with implied massive star formation 
rates of order 10 3 M Q year -1 (without correcting for possible amplification by gravitational lensing). 
The radio-to-submm spectral energy distribution for BRI 0952—0115 suggests a low-luminosity radio jet 
source driven by the AGN. 

Subject headings: radio continuum: galaxies — infrared: galaxies — galaxies: starburst — galaxies: 

evolution — quasars: individual (BRI 0952-0115, BR 1033-0327, BR 1117-1329, & 
BR 1144-0723, BR 1202-0725, LBQS 1230+1627B, BRI 1335-0417) 
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1. INTRODUCTION 



Detecting submm thermal dust emission from objects 
at z > 2 has revolu tionized our understanding of galax- 



1998 



iesat high redshift (McMahon ct al. 1994 , Hughes ct al. 
Ivison ct al. 1998a|, Smail, Iy ison, fc Blain 1997, 



Bargcr et al. 1998, Eales ct al. 199E). A number of these 



submm sources have also been detected in CO e mission 
with implied molecular gas masses > 10 10 Mq (Browr 



fc Vanden Bout 1991, Barvainis et al. 1994, Ohta et al 



1996, Omont et al. 1996a, Guilloteau et al. 1997, Frayer 



et al. 1998, 1999), and the large reservoirs of warm gas 
and dust in these systems have led to the hypothesis that 
these galaxies are starburst galaxies with massive star for- 
mation rates of order 10 3 M Q year -1 (Hughes & Dunlop 
1999, Lilly et al. 1999). 

In their study of optically selected QSOs at z > 4 from 
the APM sample (McMahon 1991, Irwin, McMahon, & 
Hazard 1991), Omont et al. (1996b) detected 240 GHz 
continuum emission in 6 out of 16 QSO's. If the detected 
emission is due to dust in the quasar host systems, the im- 
plied dust masses are « 10 s M Q . Follow-up observations 
of these dust emitting QSOs have revealed CO emission 
in three c ases so far, with i mplied molecular ga s masses ~ 
10 11 M (|Ohta et al. 1996| , |Omont et al. 1996a| Guilloteau 
et al. 1997,1999, Carilli, Menten, & Yun 1999). Omont 
et al. also point out that about half the dust emitting 
QSOs have broad absorption lines (BALs), again suggest- 
ing a gas rich environment. Omont et al. speculate that: 
"...such large amounts of dust [and gas] imply giant star- 
bursts at z > 4, at least comparable to those found in the 
most hyperluminous IRAS galaxies...". However, evidence 
for active star formation in these sources remains circum- 



stantial, primarily based on the presence of large gas reser- 
voirs, and it remains possible that the dust is heated by 
the AGN, rather than by a starburst. 

One possible method to investigate the nature of these 
high redshift dust-rich QSOs is sensitive radio continuum 
observations. A well studied phenomenon in nearby star 
forming galaxies is the radio-to-far IR correlation, i.e. the 
tight correlation found between radio continuum emission 
and thermal dust emission (Condon 1992, Helou & Bicay 
1993). The standard explanation for the radio-to- far IR 
correlation involves relativistic electrons accelerated in su- 
pernova remnant shocks, and dust heated by the interstel- 
lar radiation field. Both quantities are then functions of 
the massive star formation rate (Condon 1992, Cram et al. 
1998, Yun, Reddy, & Condon 1999, Gruppioni, Mignal, & 
Zamorani 1999), although the detailed physical processes 
giving rise to the tight correlation remain enigmatic. 

An obvious difficulty in interpreting the radio and in- 
frared emission from an AGN host system in terms of star 
formation activity is disentangling the contribution associ- 
ated with the AGN activity. A relatively minute amount of 
dust (1O 6_7 M ) heated by a luminous AGN can produce 
detectable levels of infrared and submm emission (e.g. Yun 
& Scoville 1998). Similarly, all AGNs are likely sources 
of radio emission at some level, and whether a particular 
QSO is a radio source is really a question on sensitivity. 
For example, a VLA 5 GHz radio continuum survey of 22 
optically selected high redshift (z > 3) QSOs by Schneider 
et al. (1992) has produced only one detection above the 
5cr limit of 0.2 mJy 5 while another VLA survey of Palo- 
mar Bright Quasar Survey sources at 0.02 < z < 2.1 by 
Kellermann et al. (1989) has detected 96 out of 114 (84%) 
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at a similar limiting flux level. The actual mechanism for 
generating luminosity by an AGN is not understood well 
enough nor particularly relevant to warrant further dis- 
cussions here. Given that radio luminosity of an AGN 
can span nearly ten orders of magnitude, whether the ob- 
served radio luminosity of a dusty QSO host system is 
even comparable to the value expected from the inferred 
star forming activity is our main interest. 

In this paper we present sensitive radio continuum ob- 
servations of the high redshift dust rich QSOs from the 
APM sample made with the Very Large Array (VLA) at 
1.4 GHz and 4.9 GHz. For about half the sample, our ob- 
servations are sensitive enough to detect the radio contin- 
uum associated with a starburst, as dictated by the dust 
continuum emission and the radio-to-far IR correlation. 
We find that three of the sources are consistent with the 
radio-to-far IR correlation for starburst galaxies while the 
radio emission in one of the sources is likely dominated by 
a low luminosity (FR I) radio jet driven by the AGN. We 
then examine whether other high redshift, dust rich AGN 
hosts from the literature also follow the radio-to-far IR 
correlation for starburst galaxies, and we compare these 
results to the low redshift ultraluminous IR galaxy and 
BAL QSO Mrk 231. We adopt H Q = 75 km s^ 1 Mpc" 1 
and q = 0.5. 

2. OBSERVATIONS 

The sources listed in Table 1 are from a sample of op- 
tically sele c ted z > 4 QSO s from the APM survey ( IrwirJ 
et al. 1991| , [VIcMahon 199]] , |Storric-Lombardi ct al. 1990| ) 
They arc selected for their very red color and point-like 
appearance, with the idea that the red color is due to the 
Lya forest or the Lyman break shifting into the optical 
band. An additional contribution to the red color may 
be dust obscuration within the host galaxies. Omont et 
al. (1996b) observed 16 of these sources at 240 GHz to 
sensitivity limits of 1.5 mJy or better using the bolometer 
array on IRAM 30-m telescope. They claim detections at 
> 3cr for six sources. We also include the z=2.70 QSO 
LBQS 1230+1627B, which was detected at 240 GHz by 
Omont et al. among their low redshift weak emission line 
sample. 

Observations were made with the VLA in the BnA and 
C configurations at 1.4 GHz and 4.9 GHz between 1998 
and 1999. Data from the VLA archive were also examined 
in some cases to verify our results. Integration times were 
between 1 and 2 hours long, and a total bandwidth of 2 
x 50 MHz and two circular polarizations were used. The 
absolute flux density scale was set using 3C 286 (15 Jy at 
1.4 GHz and 7.5 Jy at 4.9 GHz). After normal calibra- 
tion using VLA calibrators, the data on each source were 
edited and self-calibrated using standard methods for high 
dynamic range imaging (Perley 1999). Wide-field images 
encompassing the full primary beam were synthesized and 
deconvolved in order to mitigate image noise due to side- 
lobes from sources distributed throughout the array field 
of view. The final rms noise levels on the images and mea- 
sured source flux density are listed in columns 5 and 6 in 
Table 1. In most cases the final image noise values are 
within a factor 2 of the expected theoretical limits. The 
radio spectral index for each detected source between 1.4 
GHz and 4.9 GHz is listed in column 7 in Table 1. We 
define spectral index, a, in terms of frequency, v, and the 



observed flux density, S^, as: Sj, oc v a . 

3. RESULTS ON INDIVIDUAL SOURCES 
3.1. BRI 0952-0115 

This z = 4.43 QSO is the weakest of the claimed detec- 
tions by Omont et al. (1996b) with S 2 4 = 2.8 ± 0.6 mJy. 
Guilloteau et al. (1999) have recently detected CO(5— 4) 
emission from this QSO, confirming its gas and dust rich 
nature. The optical QSO is gravitationally lensed with a 
1" separation (Omont et al. 1996b), but the radio and mm 
observations made thus far have not had sufficient resolu- 
tion to separate the two images. 

This source is clearly detected at both 1.4 GHz and 4.9 
GHz. The source is not spatially resolved by our observa- 
tion at 4.5" x 2.7" resolution, and the measured flux den- 
sity is 460±40 AiJy and 234±27 ^Jy at 1.4 GHz and 4.9 
GHz, respectively. The radio spectral index between 1.4 
and 4.9 GHz is -0.54 ± 0.11, which is slightly flatter than 
the typical value for a starburst system, a « —0.75 ± 0.1 
(Condon 1992). 

If this were a starburst system at z — 4.43, then the 
expected 1.4 GHz flux density extrapolated from the 240 
GHz flux density would be 44 /zJy - about 5 times smaller 
than the actual observed value. Similarly, the observed 
value for af 5 ® is +0.57±0.05 while the predicted value for 
a starburst galaxy at that redshift is +0.99±0.1. There- 
fore this source has clear excess of radio emission compared 
with a normal starburst galaxy. The radio luminosity of 
this source is about 5 times that of the radio-loud ellipti- 
cal galaxy M87, and places this source in the class of low 
luminosity (Fanaroff-Riley Class I) radio galaxies. 

One way to reconcile the starburst hypothesis for driving 
the radio emission from this source is to assume that gravi- 
tational lensing preferentially amplifies the radio emission 
relative to the dust emission. It is well known that dif- 
ferential magnification by gravitational lensing can signif- 
icantly alter the apparent SED (e.g. Granato et al. 1996, 
Blain 1999a). For a pure starburst system, however, the 
radio-to-submm flux ratio should be relatively unaffected 
by lensing since radio synchrotron and ther mal dust emis- 



sion should have similar spatial extents (Condon ct al 
1991). Gravitational magnification is approximately in- 
versely proportional to angular size, and a preferential am- 
plification of radio emission is expected only if a compact 
radio nucleus is present. Several other sources discussed 
here may also be lensed, and these considerations should 
generally apply to the others as well. 

3.2. BR 1202-0725 

The z = 4.70 QSO BR 1202-0725 is the first high red- 
shift QSO where thermal dust emission was detected with 
S240 = 12.6±2.3 mJy and S350 = 49±5 mJy, respectively 
(McMahon et al. 1994). This so urce has also been d etected 



in multiple t ransitions of CO (Ohta et al. 1996 , Omont 
|et al. 1996a] , |Yun, Scoville, fc Evans 1999[ ). This source 
has the curious property that the optical QSO is a single 
source, but the mm continuum and the CO line observa- 
tions show a double source with a separation of about 4" 
dOmont et al. 19964 |Guilloteau et al. 1999]) . This dou- 
ble morphology may indicate a pair of interacting objects 
separated by only about 20 kpc. Alternatively, this may 
be a gravitationally lensed source with one of the optical 
images obscured by intervening material (see below). 
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We have detected radio continuum emission from al. 1998a , Dey et al. 1999 ) and SMM 02399—0136 ( [ivison 



BR l p02-0725 with a total flux density of 240±35 fiJy 
at 1.4 GHz. The low resolution image (Figure la; 9 = 
10.6" x 6.6") shows a slight extension of the source along 
the position angle of the mm continuum and CO sources. 
The high resolution image (Figure lb; 6 = 4.5" x 2.8") 
shows a double source with the same separation and posi- 
tion angle as the mm continuum and CO emitting sources. 
This source is also detected at 4.9 GHz with a total flux 
density of 141±15 fiJy (Figure lc). The resulting 1.4 GHz 
to 4.9 GHz spectral index is —0.43 ± 0.14, which is flatter 
than expected for a starburst galaxy. 

The observed value of a^ 5 ® for the integrated flux den- 
sity for BR 1202-0725 is +0.96±0.04, while the value pre- 
dicted for a starburst galaxy at z = 4.70 is +1.03±0.1. 
Hence the 1.4 GHz-to-350 GHz spectral index is consistent 
with a massive starburst. This conclusion is somewhat 
inconsistent with the derived radio spectral index above. 
One possible explanation for this discrepancy is free-free 
absorption of the radio continuum emission, although 
the required emission measures are large (ss fewxlO 8 pc 
cm~ 6 ). Alternatively, a flat spectrum radio AGN may be 
present, contributing significantly at 5 GHz, with steeper 
spectrum emission driven by star formation dominating 
at 1.4 GHz. This later alternative is admittedly ad hoc, 
but the situation is similar to what is seen in another 
AGN+starburst system, Mrk 231 (see below). Future ra- 
dio continuum imaging at sub-arcsecond resolution may 
help resolving this issue for BR 1202—0725. 

If the 1.4 GHz radio emission is driven by star forma- 
tion, then the required star formation rate is very large, 
« 2500 M Q year" 1 (using Eq. 1 of Carilli & Yun 1999). 
Again, given the differences among galaxies, a star for- 
mation rate derived in this way should be considered an 
order-of-magnitude estimate at best, and that the true star 
formation rate could be a factor few to ten times lower if 
the source is amplified by gravitational lensing. 

BR 1202—0725 is found in a high density environment, 
and gravitational lensing may not be needed to explain the 
double radio and submm morphology. There is a Lyman-a 
emitting companion located 2.6" northwest of the quasar 
(Hu et al. 1996), and a deep ground based K-band image 
and an HST I-band image by Hu et al. show a string of 
continuum emission extending over 4" in the same direc- 
tion. There is also a second galaxy located about 3.5" to 
the southwest of the quasar (see their Fig. 2). One of the 
two radio continuum sources is coincident with the optical 
quasar while the second radio and submm source is very 
close to but not coincident with the Lyman-a compan- 
ion. A possible faint K-band counterpart is visible in the 
image obtained by Hu et al., and the radio and submm 
source is probably also physically related to the Lyman-a 
source. Citing the close proximity (~ 20 kpc), Hu et al. 
suggested ionization of the Lyman-a companion by the 
quasar. However, external heating by the quasar can not 
explain the bright radio and submm emission in the com- 
panion. If powered internally by star formation instead, 
the star formation rate inferred from the Lyman-a lumi- 
nosity is only about IOMq yr _1 , which is about two orders 
of magnitudes smaller than inferred from the submm lu- 
minosity. Similarly weak Lyman-a emission has been seen 
in other dusty high rcdshif t submm sources such as [HR94] 
10 (ERO J164502+4626.4,|Hu k Ridgway 1994), |Cimatti et 



et al. 1998a). The 



a 



ij 4 index for the companion alone is 
T.06±0.07, consistent with the SED of a z=4.7 starburst 
system. In this scenario, the QSO is a flat spectrum radio 
source located within a massive starburst host while the 
companion is a second dust obscured starburst system. 

3.3. BRI 1335-0417 

The z = 4.40 QSO BRI 1335-0417 has been detected 
at 240 GHz by Omont et al. (1996b) with a flux density 
of 10.3±1.0 mJy and by Guilloteau et al. (1997) with a 
flux density of 5.6±1.1 mJy at 220 GHz. This difference 
in observed flux densities at these two close frequencies 
is consistent with a sharply rising submm spectral index 
(see Guilloteau et al. 1997). Given the large uncertainty 
in each of these measurements, however, we adopt a mean 
value of 8.0±2.5 mJy at 240 GHz with an uncertainty en- 
compassing both measurements. Then the expected flux 
density at 350 GHz is 31 ±9 mJy. This source has also 
been detected in CO (5-4) emission by Guilloteau et al. 
(1997), and in the CO (2-1) transition by Carilli et al. 
(1999), with an inferred molecular gas mass of » 10 11 M Q . 
Neither the CO nor the thermal dust emission is resolved 
at 1" resolution. 

The source BRI 1335—0417 is detected in the radio 
continuum at 1.4 GHz and 4.9 GHz with flux densities 
of 220±43 /iJy and 76 ± 11 fiJy, respectively. There 
is a possible confusing source about 15" to the east of 
BRI 1335-0417 (see Carilli et al. 1999). The 1.4 GHz to 
4.9 GHz spectral index is —0.85 ± 0.19, consistent with 
synchrotron emission from a starburst system. The ob- 
served value of af 4 ° for BRI 1335-0417 is +0.90±0.07 
while the value predicted for a starburst galaxy at z = 
4.40 is +1.0±0.1. Hence, the radio synchrotron and mm 
dust emission from BRI 1335—0417 are consistent with a 
massive starburst. Like BR 1202—0725, the required mas- 
sive star formation rate is very large, « 2000 M year -1 . 
Thus far there is no evidence for gravitational lensing in 
BRI 1335-0417. 

3.4. LBQS 1230+162W 

This z = 2.70 QSO is one of the low redshift QSOs in the 
Omont et al. sample with a 240 GHz detection. Omont et 
al. included this QSO in their search because they postu- 
lated that weak emission lines are indications of a possible 
dusty environment. The reported 240 GHz flux is S240 = 
7.5±1.4 mJy. The Plateau de Bure observations by Guil- 
loteau et al. (1999) have confirmed this detection, but they 
recovered only 3.3±0.5 mJy of flux at 225 GHz. The 850 
GHz flux density of 104 ± 21 mJy has also been reported 
by Benford et al. (1999), and its dusty nature appears 
secure. A previous search for radio continuum produced 
only a 3 cr upper limit of 0.24 mJy at 8.4 GHz ( Hooper et 
al. 1995|) . 



This source is detected at both 1.4 GHz and 4.9 GHz 
with measured flux density of 210±50 /^Jy and 86±11 
/zJy, respectively. The derived radio spectral index a\ 4 = 
—0.72 ± 0.22 is typical of synchrotron emission from a star 
forming galaxy. Adopting an average value of 5. Oil. 5 mJy 
at 230 GHz, the derived a? 5 4 ° is +0.82±0.07, consistent 
within the uncertainty with the predicted value for a star- 
burst at z = 2.70 of +0.71±0.10. 
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3.5. BR 1033-0327, BR 1117-1329, & BR 1144- 0723 

Along with BRI 0952-0115, these QSOs are the faintest 
sources detected by Omont et al. (1996b), with S240 
between 3.5 mJy and 5.8 mJy. The z = 4.51 QSO 
BR 1033-0327 was detected at 375 GHz by Issak et al. 
(1994) at 12±4 mJy while the extrapolated 350 GHz flux 
from the 240 GHz is 14±4 mJy. The extrapolated 350 
GHz flux for the other two sources are 16 mJy and 23 
mJy, respectively. We do not detect these three sources at 
1.4 GHz, with 4(7 upper limits between 170 /iJy and 300 
/zJy If the 240 GHz detection is real, then the 4<r lower 
limits to af 5 ® are ss +0.8 while the predicted values are 
+1.0±0.1. Hence the lower limits to af 5 4 we obtained are 
consistent with a starburst interpretation, although these 
limits by no means preclude dust heated by an AGN in 
these systems. 

4. DISCUSSION 

In the local universe, galaxies with a pronounced dust 
peak and large FIR luminosity (> 1O 12 L ) are luminous 
starburst systems fueled by a massive concentration of 
gas, enshrouded in a thick dusty cocoon (see Sanders & 
Mirabel 1996 and references therein). The 240 GHz con- 
tinuum emission detected by Omont et al. in these high 



redshift QSOs euiiespund Lu iusL fiame wavelengths of be- 
tween^ 240 pern and 350 fim, near the peak of their FIR 

dust spectral energy distribution (SED). We first consider 
their large FIR luminosity in relation to the ultraluminous 
infrared galaxies in the local universe. Then we consider 
the clear presence of optical QSOs and the possible AGN 
contribution to the radio and dust continuum emission. 

4.1. Starburst in Dusty QSOs 

The dust emission associated with these high redshift 
QSOs may be reprocessed radiation from luminous star- 
burst activity within their host galaxies is strongly sug- 
gested by the comparison of their radio to infrared SEDs 
with that of a starburst galaxy like M82 (see Figure 2). 
While there is a natural expectation of finding some ra- 
dio continuum emission associated with AGN activity in 
these QSOs, detecting radio emission exactly at the level 
expected from the dust emission has to be entirely for- 
tuitous. The fact that this trend is seen in three out of 
four QSOs detected at both radio and mm/submm wave- 
lengths is probably not coincidental. Rather, finding these 
luminous high redshift QSOs associated with gas and dust- 
rich, massive starburst systems should have been expected 
since most models of cosmic star formation history and 
galaxy formation/evolution trying to account for the faint 
submm source counts require up to two orders of magni- 
tude inc rease in comoving density of such systems a t earlier 
epochs (Guiderdoni et al. 1998 , Blain et al. 199S , Tan et 
al. 1999| ). 

The results of our radio continuum study of high red- 
shift dust emitting QSOs are summarized in Figure 3. 
This figure shows the values of af™ versus redshift for the 
seven sources discussed above. Also plotted are submm 
detected dust rich galaxies with a wide range of redshifts 
and four models for the evolution of a\^2 with redshift 
for starburst galaxies as given in Carilli & Yun (1999). 
Two of the models are semi-empirical, based on starburst 



tra of the low redshift star forming galaxies Arp 220 and 
M82. The radio-to-submm SEDs for the three sources, 
BR 1202-0725, LBQS 1230+1627B, and BRI 1335-0417, 
are consistent, within the uncertainty and the range of the 
models, with starburst-driven radio synchrotron and mm 
dust emission, despite the obvious presence of an optical 
QSO. One source, BRI 0952—0115, shows clear excess ra- 
dio emission over the expected from a starburst. 

Omont et al. have claimed > 3u detection of thermal 
dust emission at 240 GHz in 6 out of 16 QSOs from the 
APM sample at z > 4. Three of these sources have also 
been detected in CO emission, and Omont et al. point out 
that half the sources are BAL QSOs - a phenomenon that 
is typically anti-correlated with radio-loud QSOs (see Wey- 
mann 1997). Overall, the host galaxies of these QSOs ap- 
pear to be extremely gas rich, leading Omont et al. (1996a) 
to speculate that these AGN are associated with starburst 
host galaxies. However, detecting millimeter continuum 
emission alone does not serve uniquely as an indicator of 
a gas- and dust-rich galaxy vigorously forming stars since 
the presence of a flat spectrum AGN can account for the 
m m/submm detection in some cases (e.g. B 2 0902+34, 



Yun 



Scoville 1996j , |Downes et al. 1996| ). Relatively 
small amount of hot dust (10 b ~ r Af ( n) heat ed by an AGN 



can a l so produce eleva ted submm emission (Yun & Scovilk 
1998| , [Haas et al. 1998; ). The fact that six out of seven 240 



GHz detected QSOs have the same characteristic radio-to 
submm SED (including lower limits), consistent with that 
of a starburst galaxy argues that these QSOs reside in gas 
and dust rich host galaxies harboring a massive starburst. 
The gas depletion time scale is short, w 10 s years, during 
which time a significant fraction of the stars in the host 
galaxy may be formed. 

There is a considerable scatter in the data and models 
plotted in Figure 3, and this relation may only be used as 
gross redshift indicator. Understanding the nature of con- 
tributing physical causes may allow a more robust use of 
this relation. An obvious concern is the contribution from 
a radio AGN (i.e. a jet), and the presence of an ener- 
getically important AGN in general is discussed in greater 
detail below. The intrinsic scatter in the radio-to-FIR cor- 
relation is only about 0.25 in dex (see Condon 1992), much 
smaller than the scatter seen in Figure 3. The mean elec- 
tron density, which is a function of the compactness of the 
starburst region and the intensity of the starburst activity, 
can produce a varying degree of free-free absorption of ra- 
dio continuum (see Condon et al. 1991), and this largely 
accounts for the systematic difference between the empher- 
ical models based on "M82" and "Arp 220" in Figure 3. As 
discussed by Blain (1999b), dust temperature and emissiv- 
ity can also introduce variations at the submm/FIR wave- 
lengths that may contribute to the observed scatter. To 
produce a significant deviation from the ranges of models 
considered in Figure 3, however, dust temperature has to 
be exceptionally low (^ 20 K). The measured dust SEDs 
of infrared luminous galaxies in the local universe are con- 
sistently characterized by dust temperature in excess of 
30 K (see Benford 1999, Lisenfeld et al. 1999), and the 
models seem quite robust when applied to actively star 
forming galaxies. 

It is well known that the fraction of radio lou d QSOs 



galaxy relationships derived by Condon ^lyy^j, while two 
of the models are empirical, based on the observed spec- 



decreases from about 20% at z < 2 to 5% at z > 3 (McMa-. 
hon 1991] , Schneider et al. 1992j , |5chmidt et al. 1995| ). In 
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this c ase, "radio-loud" implies radio luminosity an order of 



magn: tude or more larger than the values presented herein, 
indicative of a radio jet associated with the active nucleus. 
The importance of the observations presented here is that 
the flux density limits adequate to detect the synchrotron 
radio continuum emission associated with star formation 
have been reached at least for some sources. Hence we 
are potentially probing a new population of radio-emitting 
AGN host galaxies at high redshift. 

4.2. Presence of a Luminous AGN 

There are now more than a dozen high redshift sources 
where thermal dust emission is detected. Many of these 
sources exhibiting a value of af 5 4 characteristic of a star- 
burst, and they also show clear evidence for an AGN. 
One such source is the z = 2.28 gravitationally lensed hy- 
perluminous infrared galaxy IRAS 10214+4724. Presence 
of a hidden QSO is evidenc ed by broad optical emission 
lines seen in scattered light ( Lawrence et al. 1995 , Jannuzj 
et al. 1994 ). A second source is the z = 2.903 submm 



absor ption on scales of 100 to 1 00 pc ( Bryant fc Scovilk 
1996| , [Downes fc Solomon 199S| , Parilli et al. 1998| ). The 



galaxy SMM 02399-0136, which is the first subm m iden- 
tified lumino us galaxy showing an AGN spectrum ( Ivison| 
et al. 1998a|). The third is APM 08279+5255 at z = 3.87, 



which shows an optical BAL QSO spectr um (Irwin et al 
19981 |Lewis et al. 19981 |Hines et al. 1999[) . These sources 



have also been detected in CO emission and thus harbor 



large amounts of molecular gas (Brown fc Vandcn Bout 
199lj |Frayer et al. 1998| , pownes et al. 1999)) . 

The z = 2.558 BAL QSO H 1413+117 is another submm 
and CO source (Barvainis et al. 1992,1994), but its radio 
emission has a signi ficant contribution from its radio jet 
( Kayscr et al. 199(i| ). Therefore it deviates significantly 
from the radio-submm relation for starburst galaxies as 
shown in Figure 3, similarly as BPJ 0952—0115. Other 
submm detected radio g alaxies and radio-loud QSOs such 
as 4C 41.17 at z 
8C 14354 




3.80 flDunlop et al. 1994 
635 at z 



Hughes et al 



Ivison et al 



4.25 ( Hvison 1995| . 
MG 0414+0534 at z = 2.64 (Barvainis 1998ft Tand 
MG 1019+0535 at z = 2.76 QCimatti et al. 199811 have 



comparable submm/IR luminosity, but their radio lumi- 
nosity is 2-3 orders of magnitude larger. As a result, they 
deviate so much from the radio-submm relation as to drop 
completely below the ranges plotted in Figure 3. 

While we have argued for active star formation in at 
least some of these systems, the question remains as to 
what degree the luminous AGN present in these dusty 
QSOs account for the radio and/or submm/IR luminos- 
ity. The answer to this question is difficult to obtain even 
for the ultraluminous infrared galaxies in the local universe 
(see reviews by Sanders [1999], Joseph [1999], and Genzel 
[1999] and references therein), and addressing this issue for 
the high redshift objects has to be even harder. Instead, an 
examination of a low redshift analog of a galaxy with a con- 
current AGN plus starburst in the ultraluminous infrared 
galaxy Mrk 231 may provide a useful insight. Mrk 231 is 
the most luminous infrared source at z < 0.2 (Sanders et 
al. 1988), with an infrared luminosity of 5 x 10 12 h -2 L Q . 
This is within a factor of a few to ten of the high redshift 
dust emitting sources, and Mrk 231 has been called the 
nearest BAL QSO (Forster et al. 1995). About half of its 
1.4 GHz radio emission comes from a parsec-scale radio 
jet (Ulvestad et al. 1999) while the other half is associ- 
ated with a gas disk seen in CO emission and HI 21cm 



thermal dust and non-thermal radio continuum emission 
from this disk is consistent with a massive starburst of 200 
M Q per year, and support for such a star forming disk in 
Mrk 231 comes from the tentative detection of supernova 
remnants (Taylor et al. 1999). Citing geometric and ener- 
getic reasons, Downes & Solomon (1998) have argued that 
the AGN and the starburst each account for about 1/2 of 
the total luminosity in Mrk 231. The radio flux density of 
Mrk 231 between 5 GHz and 110 GHz is a factor 3 higher 
than the starburst prediction (see Figure 2) due to the con- 
tribution from the radio jet. This component is free-free 
absorbed at lower frequencies, where the flux densities are 
more in-line with those expected for a starburst. 

Mrk 231 is an illustrative example of a dusty QSO where 
an AGN may provide significant infrared and radio lumi- 
nosity. It is the closest local analog of the dusty high red- 
shift QSOs whose radio emission has been studied here, 
but it is not certain if Mrk 231 is representative of the 
dusty QSOs in the early epochs. Nevertheless, the prop- 
erties of Mrk 231 are consistent with the sequence pro- 
posed by Sanders et al. (1988) in which merging galaxies 
evolve through a starburst phase, eventually becoming op- 
tical QSOs (see also Sanders & Mirabel 1996). The time 
scale for such evolution is of order a fewxlO 8 to 10 9 years. 
Sanders et al. propose that Mrk 231 represents a late stage 
in this process, namely the 'dust enshrouded QSO' stage. 
Based on their radio and mm properties, many dust emit- 
ting high redshift QSOs also appear to fall in this cate- 
gory although the required star formation rates are a fac- 
tor of 5 to 10 larger than IR luminous galaxies at low 
redshift (in the absence of amplification by gravitational 
lensing). Clearly, deeper radio and mm observations at 
sub-arcsecond spatial resolution are required to test this 
hypothesis. Such observations will become possible with 
the upgraded VLA at cm wavelengths and with the future 
Atacama Large Millimeter Array in Chile. 

Despite the large error bars associated with the data 
points and a significant spread among the models, the 
dusty AGNs (filled and empty circles) as a group may 
suggest a systematic deviation by about 0.2-0.3 from the 
radio-submm relation obeyed by the pure starburst sys- 
tems (empty squares) in Figure 3. It is tempting to inter- 
pret this trend as a residual effect of energy input by an 
energetic AGN in some cases, similarly as in Mrk 231 and 
to a lesser degree as in H 1413+117 and BR 0952-0115. 
Better and additional data on the high redshift dust-rich 
sources identified by SCUBA as well as more refined mod- 
els with improved understanding of the intrinsic scatter 
are needed to examine the validity of this possible trend 
in detail. 

5. SUMMARY 

The 1.4 and 4.9 GHz radio continuum emission is de- 
tected in four out of seven high redshift QSOs that had 
previously be detected in thermal dust emission at 240 
GHz by Omont et al. (1996b). In three of the cases, the 
cm and mm emission are consistent with starburst activity. 
One of the QSOs clearly hosts a radio AGN. If the detected 
radio continuum is interpreted solely in terms of AGN ac- 
tivity instead, the observed luminosity being within a fac- 
tor of 2 of the value inferred for starburst activity would 
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be fortuitous, although such a possibility cannot be ruled 
out. 
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Caltcch, under contract with NASA. The National Ra- 
dio Astronomy Observatory (NRAO) is a facility of the 
National Science Foundation, operated under cooperative 
agreement by Associated Universities, Inc. 
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Fig. 1.— (a) A low resolution 1.4 GHz image of the z = 4.70 QSO BR 1202-0725 made with the VLA. The angular resolution is 10.6"x6.6" 
(FWHM, with major axis position angle = 90°). The crosses show the positions of the mm continuum and CO sources by Omont et al. 
(1996b). The contour levels are: -70, 70, 105, 140, 175, 210, and 245 /iJy beam" 1 , (b) The same as la, but now at higher resolution: 
4.5" X 2.8" (FWHM) with major axis position angle = 82°. The contour levels are -64, -32, 32, 64, 96, and 128 /iJy beam -1 . The emission 
appears to be resolved into two separate sources roughly coincident with the positions of the two submm sources reported by Omont et al. 
The flux densities and J2000 positions of the two sources arc: 70±35 /iJy at 12 h 05 m 23.16 s , -07° 42' 32.4" and 133±35 /iJy at 12 h 05 m 
23.02 s , -07° 42' 29.8" with position uncertainties of ± 1". (c) A 4.9 GHz image of the z = 4.70 QSO BR 1202-0725 made with the VLA. 
The angular resolution is 16.2" X 7.2" (FWHM, with major axis position angle = -69°). The contour levels are: -45, -30, 30, 45, 60, 75, 90, 
105, 120, and 135 ^Jy beam" 1 . 

Fig. 2. — The spectral energy distribution of 1335—0415 (squares) and BR 1202—0725 (circles) are compared with that of the prototypical 
starburst galaxy M82. The vertical scale is set for the flux density of M82, and the SEDs for all other sources are normalized to M82 at the 
peak of their dust SED. The intrinsic 1.4 GHz luminosity (i.e. in the source rest frame) of BRI 1335—0415 and BR 1202—0725 arc a factor 
1500 or so higher than M82. The SED for Mrk 231, an ultraluminous infrared galaxy and the nearest BAL QSO, is shown (stars with an 
error bar) as an illustrative example of a dusty QSO hosting a radio AGN. It has about 100 times larger 1.4 GHz luminosity than M82, and 
its flat spectrum radio nucleus and a radio jet accounts for a systematic offset between 10 10— 11 Hz. 

Fig. 3. — The data points show derived values of af 5 ® versus redshift for the sources in Table 1, including 4<r lower limits. The lines are 
models derived for starburst galaxies from Carilli & Yun (1999). Two power-law models for a^2 ^ or star f° rmm g galaxies as a function of 
redshift are shown with a solid curve (a su bmm = +3.0) and a short dashed curve (a su bmm = +3.5). The assumed radio spectral index is 
—0.8. The dotted and long-dashed curves show the expected values of aJ50 versus re dshift based on the observed SED of the star forming 
galaxies M82 and Arp 220, respectively. The Omont et al. sample QSOs and Mrk 231 arc plotted as filled circles while the lower limits for the 
three undetected Omont et al. QSOs are shown as arrows. Dusty sources from the literature with clear evidence for an AGN are also plotted 
as empty circles while the sources without evidence for an AGN are plotted as empty squares. 
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Table 1. Measured Fluxes and Radio-to-Submm Spectral Index 



Source 


z 


RA (J2000) 


Dec (J2000) 


SlAGHz 

(^Jy) 


SsGHz 

(/^Jy) 


a\ 4 (obs) 


(^Jy) 


af 5 4 °(obs) b 


af 5 4 °(pred) c 


BRI 0952-0115 


4.43 


09 ft 55 m 00M2 


-01°30'06".9 


460 ± 40 


234 ± 27 


-0.54 ± 0.11 


11±3 


0.57±0.05 


0.99±0.10 


BR 1033-0327 


4.50 






<170 






14 ±4 


> 0.79 


1.00±0.10 


BR 1117-1329 


3.96 






<170 






16±5 


> 0.83 


0.93±0.10 


BR 1144-0723 


4.14 






< 300 






23±6 


> 0.78 


0.95±0.10 


BR 1202-0725 


4.70 


12 ft 05 m 23 s .09 


-07°42'31".l 


240 ± 35 


141 ± 15 


-0.43 ±0.14 


50±8 


0.96±0.04 


1.03±0.10 


LBQS 1230+1627B 


2.70 


12 ft 33 m 10 s .32 


+ 16°10'50".8 


210 ± 50 


86 ± 11 


-0.72 ±0.22 


20±7 


0.82±0.07 


0.71±0.10 


BRI 1335-0417 


4.40 


13 ft 38 m 03 s .44 


-04°32'34".l 


220 ±43 


76 ± 11 


-0.85 ±0.19 


31±9 


0.90±0.07 


1.00±0.10 



a Source flux density at 350 GHz (850 //m) is extrapolated from the 240 GHz flux density by Omont et al. (1996b) assuming a spectral 
index of ±3.25. For BR 1033—0327 and BR 1202—0725, the 800 //m measurements by Isaak et al. (1994) are used instead. 
b Observed 1.4 GHz to 350 GHz spectral index for each source or 4<r lower limits. 

c Predicted value for a starburst galaxy using Eq. 3 of Carilli & Yun (1999). The uncertainty of 0.1 reflects a wide range of values 
predicted by the different models. 



References. — (1) Omont et al. (1996b); (2) Issak et al. (1994); (3) McMahon et al. (1994); (4) Guilloteau et al. (1997); (5) Guilloteau 
et al. (1999) 



